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Abstract. Drilling holes in a bulk high-Tc superconductor enhances the oxygen 
annealing and the heat exchange with the cooling liquid. However, drilling holes also 
reduces the amount of magnetic flux that can be trapped in the sample. In this paper, 
we use the Bean model to study the magnetization and the current line distribution 
in drilled samples, as a function of the hole positions. A single hole perturbs the 
critical current flow over an extended region that is bounded by a discontinuity line, 
where the direction of the current density changes abruptly. We demonstrate that the 
trapped magnetic flux is maximized if the center of each hole is positioned on one of the 
discontinuity lines produced by the neighbouring holes. For a cylindrical sample, we 
construct a polar triangular hole pattern that exploits this principle; in such a lattice, 
the trapped field is ^ 20% higher than in a squared lattice, for which the holes do not 
lie on discontinuity lines. This result indicates that one can simultaneously enhance 
the oxygen annealing, the heat transfer, and maximize the trapped field. 
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1. Introduction 

High-temperature bulk superconductors are very promising materials for permanent 
magnet applications [H [21 [Sj HI [5]. They can be used in magnetic bearings (in the 
Maglev train [6j or in frictionless hnear translation systems [7j) and in rotating machines 
(synchronous motors [8l[9] or flywheels for energy storage p^fTT] ). At the liquid nitrogen 
temperature, such magnets are able to trap up to 3 T [T2]. When cooled down to 29 K, 
the maximum trapped field can reach 17 T [T3] . 

Recently, it has been proposed to drill arrays of columnar holes inside high-Tc 
superconducting magnets in order to improve their chemical and thermal properties [TU 
[T5] . First, the holes reduce the oxygen diffusion wall and enhance the oxygen annealing 
process [16j. Second, the larger exchange surface increases the heat transfer with 
the environment and is thus beneficial for the cooling of the superconductor [171 . A 
rapid cooling is required for instance when a superconductor is magnetized with a 
pulsed field p^, because the dissipative motion of vortices tends to raise rapidly the 
temperature of the material and thus to reduce both the critical current density and 
the trapped magnetic flux. A third (although counterintuitive) advantage to drilling 
holes in a superconductor is to improve their mechanical properties. Samples can be 
strengthened by impregnating the holes with a reinforcement resin that prevents cracks 
from developing [19], for instance as a result of strains induced by the Lorentz force |13j . 

Drilling holes in a superconductor is however detrimental to its magnetic properties. 
It was found in [201 m\ that removing superconducting matter decreases both the full 
penetration field and the trapped flux. Holes also lead to macroscopic changes in the 
current distribution. In the Bean model, the current stream lines near a hole abruptly 
change their direction along discontinuity lines [23] and circle the hole in a region 
that extends far beyond the hole itself. This effect is enhanced in thin films, as the 
magnetic flux density displays sharp peaks at the discontinuity lines. Such macroscopic 
changes of the magnetic flux were observed with magneto-optical imaging of thin films 
with macroscopic defects [211 [2S]- For bulk samples, studies based on the Bean model 
already pointed to the magnetization drop that results from drilling holes [20] ■ It was 
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also shown that for a given lattice, the magnetization drop increases with the diameter 
of the holes [20l|21]. It has been measured in [21] that increasing the hole diameter by a 
factor of 2 results in a magnetization drop of ~ 80%. In the particular limit of YBCO 
thin films of rectangular shape with microscopic holes, the Bean critical state has also 
been simulated in [22]. However, to our knowledge, none of these previous works has 
studied the influence of the hole pattern on the magnetization drop. 

In this paper, we investigate the effect of the arrangement of holes on the 
magnetization drop of drilled samples, by studying the current distribution and the 
interaction among the influence regions of the holes. For that purpose, we develop 
an algorithm based on the Bean model and on an observation made by Campbell and 
Evetts [23] to calculate the magnetic field in the critical state for an infinitely long 
drilled sample with an arbitrary cross section. 

This paper is organized as follows. The algorithm is discussed in section [2] and is 
used in section |3] to calculate the magnetization of a sample with a semi-infinite cross 
section and a single hole. In sectionHl we study the magnetization of samples with either 
two or three holes, as a function of their relative positions. Section [S] is devoted to the 
magnetization drop in samples with either a semi-infinite or a circular cross section and 
holes disposed on a lattice pattern. Section [H] concludes this work. 

2. Model for the magnetic field distribution in drilled samples 

In this paper, we neglect demagnetization effects and focus on superconducting samples 
that are infinitely long and have either a semi-infinite or a circular cross section. Figure 
[1] shows a sample with a semi-infinite cross section containing a circular hole of radius 
R located at a distance D from the border. The cross section lies in the x — y plane. 
We further assume that the applied magnetic field Ha is oriented along the 2;-axis and 
is uniform. Its amplitude is such that Hd ^ Ha ^ Hc2- We assume strong pinning 
and neglect surface barrier effects. Under these assumptions, the distribution of the 
magnetic field in the sample cross section is described by the Bean model ^6], which 
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Figure 1. Sketch of an infinitely long sample with a semi- infinite cross section drilled 
by a single hole of radius R located at a distance D from the border. The flux front 
can reach the point P by following two paths: directly from the border, with a path 
length c?i, or via the hole, which acts as a radial source, yielding a path length d2- 
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where Jc is constant, while i represents the distance traveled by the flux front to reach 
a given point, P. 

Let us illustrate the procedure to determine the field distribution in the example 
shown in Figure [H For a point P located at a distance di > D, the flux front can travel 
along two paths: it can reach P directly from the border, with a path length di, or via 
the hole which acts as a radial source of magnetic field [23], with a total path length 
d2- Following Campbell and Evetts ^23j, we assume that the flux front travels along the 
shortest path. Hence, the points in a flux front are located at a fixed length i = Ha/ Jc 
from the border, where i is evaluated as di or d2, whichever is smaller. The magnetic 
field, H = Hz, at a given point, P, can then be calculated by determining the length, 
i, of the shortest path that reaches P and by evaluating 

H = Ha-Jct (2) 

Once the distribution of H is known, the current stream lines can also be easily obtained, 
as they coincide with the lines of constant magnetic field. Finally, the magnetization of 
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Figure 2. Simulated magnetization curve of an infinitely long sample with a circular 
cross section (radius a) containing one hole of radius R — 0.05 a located at a distance 
D = 0.2 a from the border. The applied field, Ha, first increases from to 2H* , where 
H* is the penetration field, then decreases from 2H* to —211*, and finally increases 
again to H*. The magnetization obtained for Ha > H* is equal in magnitude to the 
remnant magnetization. 



where S is the sample cross section (for a semi-infinite cross section, S is arbitrarily 
chosen to be a square section of unit length) and we assume B = HqH. 

We use these principles to study samples with an arbitrary number of holes and 
construct an algorithm that calculates the magnetic field distribution as a function of the 
hole radii and positions. In the following sections, we address two questions: (i) what 
is the magnetization of a given sample that is subjected in the zero field cooled state 
to an applied field. Ha-, and (ii), what is the remnant magnetization that is obtained 
when the same sample is first magnetized above twice its penetration field before the 
applied field returns to zero. For samples with a finite cross section, both magnetizations 
actually have the same magnitude (they have opposite signs, however), provided that 
the applied field in situation (i) is larger than the penetration field. Figure [2] illustrates 
this equivalence for the case of a sample with a circular cross section containing one 
hole, a case which will be treated in section [51 Such equivalence cannot be found for 
samples with a semi-infinite cross section, as they are never fully penetrated. We will 
nevertheless consider these systems when subjected to an increasing field (case (i)), 
because these situations allow us to understand the interaction between different holes. 



the sample is given by 




(3) 
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3. Samples with one hole 

3.1. Current lines 

We first consider a sample with a semi-infinite cross section drilled by a single hole of 
radius R located at a distance D from the border, and apply a magnetic field Ha in the 
zero field cooled state. Following the main principles of our algorithm, we know that 
the magnetic field can reach a given point by two distinct penetration routes. We can 
thus identify two regions: one for which the direct penetration from the border has the 
shortest path, and one for which the radial penetration via the hole has the shortest 
path. Hence, the boundary between these regions is characterized by the equality of 
path lengths, 

X = D + r - 2R (4) 

where x is a cartesian coordinate along an axis that is perpendicular to the external 
boundary and r is the distance from the hole center to the point where we determine 
the path lengths (see Figure [3]). We thus find that the boundary defined in (jlj) is the 




Figure 3. Simulated current lines (or constant magnetic field lines) in a sample with 
a semi-infinite cross section and a single hole. The hole has a radius R — 0.05 and it is 
located at distance D — 0.2 from the border. Here, the unit length corresponds to the 
length of one of the segments of the dashed contour. The applied magnetic field is such 
that the penetration depth is given as ^dcpth — 0.6. For a point Q on the discontinuity 
line, the path lengths of both penetration routes, x and D — 2R + r, are equal. 
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locus of points for which the difference between the distance to the external boundary 
and that to the hole center is equal to a constant, D — 2R. This locus is a parabola 
whose vertex is located at {x, y) = {O, D — R), whose directrix runs along y = D — 2R, 
and whose focus lies at {0,D). In cartesian cooordinates, the parabola equation reads 

x = ^ + D-R. (5) 

It is plotted as a thick line in Figure [3] for the case of a hole of radius R = 0.05 located 
at a distance D = 0.2 from the border. Here, the unit length corresponds to the length 
of one of the sides of the square delimited by the dashed contour; all the distances are 
normalized to this length. 

Equation also characterizes the current discontinuity line. As explained in the 
previous section, the current stream lines can be constructed from the contour lines 
of constant magnetic field. These lines follow straight segments outside the parabola, 
where the distance to the border is the shortest, and arcs of circle inside the parabola, 
where the penetration path through the hole is the shortest. The current lines abruptly 
change their direction on the parabola, which is thus a discontinuity line. Figure [3] 
shows the current lines obtained when the field is applied in the zero field cooled state 
and is raised to a finite Ha- In the particular case shown, the applied field corresponds 
to a penetration length £depth = 0.6. 

3.2. Influence of the hole radius on the magnetization drop 

In the case of a sample containing one single hole, the magnetization can be calculated 
in two ways: either numerically, by using the algorithm described in section [21 or 
analytically, by calculating the magnetic flux inside and outside the parabola of 
Equation ([5]). The relative magnetization drop incurred by the drilled sample is then 
given by 

AM ^ |Mo-Mi| 

|Mo| |Mo| ' ^ ^ 

where Mq is the magnetization of a sample without a hole and Mi is that for a sample 
with a single hole. The calculations are carried over a square of unit side length. In the 
particular case considered, the hole center is located at a; = 0.2 and y = 0, and we let 
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Figure 4. Magnetization drop in a sample with one hole and a semi-infinite cross 
section as a function of the hole radius. The hole is located at a distance D=0.2 
from the border. The sample is limited to a unit surface for the calculation of the 
magnetization. Here, Mo is the magnetization of a sample with the same geometry 
and no hole. The applied magnetic field is Ha = 0.6 Jc- The solid line corresponds to 
the analytical solution, and the squared symbols to that obtained with the numerical 
algorithm. 



the magnetic field penetrate up to a length £depth = 0.6. In units of Jc, the applied field 
is thus given as Ha = Jc^dcpth = 0.6 Jc- These choices guarantee that the flux does not 
extend further than x = 1 in the hole influence region. 

Figure m shows the relative magnetization drop, AM/|Mo|, as a function of the hole 
radius R (the analytical result is shown with solid lines, the numerical one is plotted 



with square symbols). The detailed analytical calculations are given in Appendix A 



Analytical and numerical calculations are in good agreement. We observe that AM/\Mq\ 
increases with the radius of the hole, as expected intuitively and illustrated in Hall probe 



mapping experiments [20] . As shown in Appendix A , a series expansion of the analytical 
result for the magnetization drop around R = yields 
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This is not a trivial result! One could have naively expected that the magnetization 
drop roughly scales either as the area of the hole, AM oc R^ -, or as the area of the 
region delimited by the parabola, AM oc \fR- From Equation ([7]), we conclude that an 
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intermediate situation occurs. 



4. Samples with several holes 
4.I. Samples with two holes 

Consider now a sample with a semi-infinite cross section and two holes. Both holes have 
the same radius R = 0.05 and are separated by a constant distance d = 0.2. The first 
hole is placed at a distance D = 0.2 from the border. Again, the semi-infinite surface is 
limited to a square whose sides have a unit length. Let us vary the angular separation 
6 between the holes and study the interactions between the regions of influence of the 
holes; their interaction should depend on whether the center of the second hole lies inside 
or outside the parabolic influence region of the first hole. We work in polar coordinate 
and take the origin at the center of the first hole. The coordinates of the second hole 
are given by (r, 6) with r = d. The discontinuity line produced by the first hole has the 
equation 



Thus, the center of the second hole lies on the parabola when r = 0.2, and hence when 
e = 30°. 

Figure [5] shows the current lines for four different angular positions 6. For 6 = 10°, 
the center of the second hole is located outside the influence region of the first hole. 
A new discontinuity parabola appears around the second hole. The two parabolic 
curves merge between the holes and form a common discontinuity line. This last line 
corresponds to the locus of points for which the difference between the distances to 
each hole center is equal to a constant; the discontinuity line is therefore a branch of a 
hyperbola. When 6 increases further, the second hole is pushed away from the border 
and, for 6 = 30°, enters the influence region of the first hole. Again, each hole produces 
a parabolic discontinuity line and the two lines merge into a branch of hyperbola. As 
the second hole goes deeper in the region of influence of the first one, the hyperbola 
opens up. The surface of the combined region of influence of the holes increases with 6 
and reaches a maximum for 6 = 90°. 



Drilled HTS: how to arrange the holes to maximize the trapped magnetic flux ? 



10 




Figure 5. Simulation of the current lines in a sample with a semi- infinite cross section 
and with two holes. The holes have a radius R = 0.05. The unit length corresponds to 
the side of the square delimited by the dashed lines. The first hole is located at D = 0.2 
from the border. The separation distance between the holes is constant, d = 0.2. The 
second hole is located at = 10° {&), 6 ^ 30° (b), 6* = 60° (c) and 6 = 90° (d). The 
thick lines represent the discontinuity lines. 

We evaluated the magnetization drop induced by the second hole as 
AM _ IM1-M2I 

|Mi| |Mi| ' ^ ' 

where Mi is the magnetization for the sample with hole 1 only, and M2 is that for 
the sample with holes 1 and 2. We can in principle evaluate this expression either 
by following the numerical method exposed in section [21 or analytically. However, 
analytical calculations rapidly become tedious when several holes are involved; we will 
thus restrict ourselves to numerical results from now on. The magnetization drop is 
plotted in Figure [6] as a function of the relative angular position of the holes. For 
small 6, the magnetization drop decays as the angle is increased. This result follows 
from the fact that the second hole is pushed away from the border as 6 increases; the 
hole is thus threaded by a lower magnetic flux and its effect is reduced. By contrast, 
for large angles, AM/|Mi| increases with 6 because the influence region of the second 
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Figure 6. Magnetization drop in a sample with two holes and a semi-infinite cross 
section, as a function of the relative angular position of the holes. The first hole 
has a radius R = 0.05 and is located at a distance D = 0.2 from the border. Here, 
the reference magnetization is calculated on a sample with the same geometry, but 
containing only the first hole. The second hole has a radius R ~ 0.05 and is located at 
a constant distance from the first hole {d ~ 0.2). The applied magnetic field is given as 
Ha = 0.6 Jc- The angular positions (a), (b), (c) and (d) correspond to the situations 
in Figure [5l We observe that AM/|Mi| has a minimum when the center of the second 
hole lies on the discontinuity line of the first hole. 

hole becomes larger. Hence, the optimal position of the center of the second hole is 
right on the discontinuity line of the first hole: the magnetization drop is minimum 
for 6 = 6opt = 30°. Although the results are not shown, we have also studied the 
situation with a fixed relative angular position and a variable separation distance and 
found similar conclusions. 

Some insight on these results can be gained by examining how the flux penetrates 
the system. The flux front near the second hole is sketched in Figure[7]for 6 = ^opt = 30° 
(a), ^ < ^opt (b), and 6 > Oopt (c). We can observe that the flux front reaches the second 
hole tangentially in all cases. However, for 6 = 6'opt, the flux front is tangent to the 
second hole simultaneously in the region inside the discontinuity parabola produced by 
the first hole (circular front) and in the region outside the parabola (straight front). 
The simultaneous penetration from the two regions appear to be necessary for reducing 
the effect of the second hole on the magnetization of the sample. 
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Figure 7. Sketch of the flux front tangent to hole 2 for angular positions 6 ~ 6'opt 
(a,), 9 < Oopt (b), and 6 > dopt (c). The dashed parts are the remnant flux front in the 
influence region of the second hole. When the center of the second hole is located on 
the discontinuity parabola, the flux front is tangent to the second hole simultaneously 
in the regions inside and outside the parabola. 



4-2. Samples with three holes 

We now turn to adding a third hole to the optimized two-hole pattern of Figure [7]- (a), 
where the first hole is again located a.t D = 0.2 away from the external border and 
the center of the second hole lies on the discontinuity parabola of the first one, at a 
distance d = 0.2. The three holes have a radius R = 0.05. The third hole is located at a 
constant distance d = 0.2 from the second hole. Adding a third hole reduces again the 
magnetization of the sample. The magnetization drop is now given by 



where M2 is the magnetization of a sample containing only the first two holes, and M3 
is the magnetization for a sample containing three holes. 

Figure [H] shows the magnetization drop, AM/IM2I, as a function of the angular 
position of the third hole; the inset shows the configuration of the holes. The center of 
the third hole can be located in three different regions: region A, a region that is not 
affected by holes 1 and 2, and regions B and C, that respectively correspond to the 
influence regions of hole 2 and hole 1. Consider first that hole 3 lies in region A. As 
6 increases, the magnetization drop, AM/IM2I, decays because the distance of hole 3 
from the border decreases (such a behaviour was already observed with the two-holes 
pattern). The magnetization drop is minimum when the center of hole 3 lies on the 
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Figure 8. Magnetization drop of a sample with three holes and a semi-infinite cross 
section, as a function of the relative angular position between the second and the third 
hole. The first hole has a radius R = 0.05 and is located at a distance D = 0.2 from 
the border. The second hole has the same radius and is located on the discontinuity 
parabola of the first hole, at a distance d = 0.2. The reference magnetization is 
calculated on the same sample with only the first and the second holes. The third hole 
has a radius R = 0.05 and is located at a constant distance d = 0.2 from the second 
hole. The applied magnetic field is Ha = 0.6 Jc- Again, the maximum magnetization 
is obtained when the center of the third hole lies on the discontinuity parabola of the 
second one. 

discontinuity parabola of the hole 2, separating regions A and B. U 9 increases further, 
hole 3 enters region B, and AM/IM2I increases again, to reach a maximum when the 
center of hole 3 reaches the remnant parabola of hole 1 (this hne does not appear as a 
discontinuity line in the current line distribution). Then, AM/IM2I decreases until hole 
3 reaches the boundary between regions B and C (this discontinuity line is a hyperbola). 
As it continues through region C, AM/IM2I increases again. We can thus conclude that 
the magnetization drop is minimized each time the center of the hole is located on a 
discontinuity line. Note however that the values of AM/\M2\ on a minimum arc not 
equal; the lowest value of AM/IM2I is achieved on the boundary between regions A and 
B. 
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5. Influence of the type of lattices 

5.1. Sample with a semi-infinite cross section 

Consider applying the results of the previous sections to construct a lattice containing 
many holes. Let us first compare the magnetization for two lattices. The first lattice is 
a squared lattice, where a line of equidistant holes (separated by a distance d) parallel 
to the external boundary is reproduced periodically every distance d in the direction 
perpendicular to the border. The second lattice is obtained from the squared lattice 
by shifting every other row by half the length of a unit cell, leading to a particular 
realization of a centered rectangular lattice. The two lattices are represented in the 
insets of Figure [91 Since we are working with semi-infinite cross sections, we work in 
a square of unit length, which we arbitrarily choose to contain five rows of holes. The 
holes have a common radius, R = 0.05, and the lattice constant is fixed aX d = 0.2. 
Note that the hole density is equal for the two lattices. 




cm D.tH fl.OI 003 D.CM 0.05 0.06 D.07 0.09 
RadluiR 

Figure 9. Magnetization drop (as compared to a sample with the same geometry and 
without lioles), as a function of the hole radius in a sample with a semi-infinite cross- 
section (limited to a square of unit length) with two different lattice configurations. 
The open circle symbols correspond to a squared lattice with a lattice constant d — 0.2 
and the filled squared symbols refer to a centered rectangular lattice with the same 
lattice constant. The insets show the lattices. The number of holes for a given radius 
is the same in each lattice. The applied magnetic field is Ha — 0.6 Jc- The centered 
rectangular lattice produces the largest magnetization. 
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The magnetization drop is defined as 
AM _ |Mo - M| 

|Mo| " |Mo| ' ^''^ 
where M is the magnetization of the sample and Mq stands for the magnetization for 

a sample with the same geometry but without holes. The applied field is carefully 

chosen to be Ha = 0.6 Jc so that the flux front stays within the square of unit length. 

AM/ 1 Mo I is plotted as a function of the hole radius in Figure [H We find that the 

centered rectangular lattice produces a larger magnetization than the squared lattice. 

Although not shown, we checked that this result is independent of the hole separation 

distance d. This result naturally follows from the conclusions of the previous sections: 

in the centered rectangular lattice, the holes are located on the discontinuity parabola 

of the neighbouring holes and the magnetization is maximized. 

5.2. Sample with a circular cross section 

Consider next infinitely long samples with a circular cross section. These samples have 
a geometry which is more realistic for bulk HTS applications. The Bean model in 
infinitely long geometries describes well the magnetic properties in the median plane of 
a cylinder with a finite height, provided its height h is large with respect to its diameter 

D PZIEH]. 

We found earlier that in a centered rectangular lattice, the holes were placed on 
the discontinuity lines of the neighbouring holes. This placement helped increasing 
the magnetization. However, this result is no longer correct for circular cross sections, 
because of the flux front geometry. The flux front is now circular, and as flux penetrates 
the system, the critical currents flow around concentric circular trajectories. Such a 
geometry is not compatible with the symmetry imposed by a centered rectangular 
lattice. The current lines for a centered rectangular hole pattern are represented in 
Figure [TUl- (a). One can observe for instance that the hole indicated by the arrow is not 
located on a discontinuity line. 

We can construct another lattice, that uses the circular shape of current lines and 
places the holes on discontinuity lines. Figure [TU]-(b) shows such a realization, which 
we name a "polar triangular lattice". The holes are positioned on concentric layers 
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separated by a distance d. Inside each layer, the hole have a common angular separation. 
Furthermore, the holes are shifted every other layer by half their angular separation. 
This ensures that the holes are located on discontinuity lines. 

By comparison with the squared lattice, one could also define a polar squared 
lattice where the hole angular position is not shifted from one layer to the next. The 
polar squared lattice, the polar triangular lattice, the squared lattice and the centered 
rectangular lattice are represented in a sample with a circular cross section in Figure fTO]- 
(c). Each lattice contains sixty holes with a radius R = 0.05 a, where a is the radius of 
the cylinder. The lattice constant is fixed to d = 0.2 a and, for the polar lattices, the 
angular separation within a layer is fixed to ^ = 20°. The corresponding magnetization 
drops are shown in Figure [TTJ The applied field is such that the cylinder is fully 
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Figure 10. (a) Simulation of the current lines in a sample with a circular cross section 
and a centered rectangular hole lattice. The arrow indicates a hole which is not located 
on a discontinuity line, (b) Simulation of the current lines in a sample with a circular 
cross section and a polar triangular hole lattice, (c) Lattice pattern under consideration 
in samples with a circular cross section. 
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Figure 11. Magnetization drop in a sample with a circular cross section (unit radius) 
and sixty holes arranged in the four lattices presented in Figure [TDKc). The holes have 
a radius R — 0.05 a, the lattice constant is fixed at c? = 0.2 a, where a is the radius of 
the cylinder. The angular separation in the polar lattice types is 20°. The sample is 
fully penetrated, Ha = Jc a- The sample with the polar triangular hole lattice shows 
the smallest magnetization drop. 



penetrated, = a. The reference magnetization Mq is calculated for a sample 
with the same geometry and without holes. We thus find that the sample with the 
polar triangular hole lattice, which aligns holes of each layer on the discontinuity lines 
produced by the previous layers, has the smallest magnetization drop. According to the 
arguments of section [21 this lattice will also have the highest trapped field. 

The results of this study are based on the neglect of demagnetization effects and 
on the assumption that the critical current density is independent of the magnetic field 
strength. However, it is worth mentioning that, under the hypothesis of a constant 
critical current density, the remnant magnetization per unit volume is not influenced by 
demagnetization effects. Therefore, the result produced by the Bean model is also valid 
for a cylinder with a finite height, as already observed for bulk cylinders in Reference |29j . 
Thus, the conclusions drawn about the maximum magnetic flux that can be trapped 
remain applicable for cylinders of finite height. 

6. Conclusions 



The magnetization drop induced by the removal of superconducting material in drilled 
samples has been studied numerically for different hole arrangements. We have 
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developed an algorithm which calculates the magnetization in the critical state of 
infinitely long samples with an arbitrary hole pattern. The main principle of this 
algorithm lies on the shortest travel path for the flux front to reach a given point in the 
cross section. The algorithm successfully reproduces the discontinuity parabola attached 
to a single hole. The dependence of the hole radius on the magnetization drop of a sample 
with one hole indicates that the loss in magnetization scales neither with the surface of 
the hole (oc R'^), nor with the surface of its parabolic region of influence (oc \/R), but as 
a surface of intermediate size, that is as li^^"^. From the simulations of samples with two 
and three holes, we have shown that in order to maximize the magnetization, the holes 
should always be located on discontinuity lines of their neighbours. The optimal lattice 
arrangement aligns the holes on the discontinuity lines and depends on the sample cross- 
section: we obtained the largest magnetization with a triangular hole lattice for samples 
with a semi-inflnite cross section, and with a polar triangular lattice for samples with a 
circular cross section. 
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Appendix A. Analytical calculation of the magnetization drop in a sample 
with a semi-infinite cross section and one hole 

Let us consider an infinitely long sample with a semi-infinite cross section and one hole 
of radius i?, located at a distance D from the border. The current line distribution 
is given by Figure lAH -(a). The magnetization drop is defined as in the main text as 
AM/ 1 Mo I = |Mo — Ml I /I Mo I, where Mo and Mi respectively are the magnetization 
of the sample of unit surface, with and without a hole. The applied magnetic field 
oriented along the z-axis is given by Ha and the penetration depth is £ = ^depth- The 
magnetic field has only a single component oriented along the z-axis. The magnetization 
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Hint 
bnglli 





Figure Al. (a) Sketch of the current line in a sample with a semi- infinite cross 
section drilled by a hole of radius R and located at a distance D from the border, 
(b) Geometrical representation of the surface where we calculate magnetic flux for the 
magnetization difference AM. 



is calculated as 
M 



^ [ BdS-Ha = —B- Ha 

Lno Js /in 



(A.l) 



fioS Js /io 
with S* = 1 is the cross section of a square of unit length. 

As the presence of the hole only modifies the flux front and the current lines inside 

the parabolic discontinuity line, the magnetization difference AM can be decomposed 

as 

Mo - Ml = — fi3" - + B^)) (A.2) 

where B"'''^''^ are respectively the average magnetic flux evaluated in the grey areas 
represented in Figure IXTI -(b). The three contributions are 

S'^ = ^/ / ^ B{x)dydx, (A.3) 

O JD-R J-2VRx-Dx+R^ 
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where B{x) = /io(4epth - x), 

S' = l r /""/^olWh -{D- R))rdrde, (A.4) 
o Jo Jo 

""^^ 7depth-p-i?)), (A.5) 



= - <! / / B{r)rdrdd + / B{r)rdrd9 \ , (A.6) 



s 

and 

2 

(S I Jr Jdmax JR J 

where B{r) = /io(^depth — r — D + 2R). The integrals in the right side are carried in 
polar coordinates, with the origin fixed at the center of the hole. In these coordinates, 
the parabola equation is given as 

- 

^parabola Z 77 • 

1 — COS 

The fiux front intersects the parabola at an angle Omax that is defined by 

2R 



(A.7) 



1 — cos t/ri 

with 



(A.8) 



Tmax = Wh -D + 2R. (A.9) 

Carying the angular integrals flA.6p . we arrive at 

-B'^ = ^ ^™a.(4cpth -r-D + 2R)rdr (A.IO) 

+/io l^"""""^ (vr - e^^ax) (Wh -r-D + 2R)rdr^ (A.ll) 

The integral over r is computed numerically. 

The series expansion around i? = of AM = \Mq — Mi\ yields 
Jc /32 



AM = ^-^R^(i^^^,^ - Df + 0(i?5/2) j (A_^2) 
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